Colloids dragged through a polymer solution: experiment, theory and simulation 
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We present micro-rheological measurements of the drag force on colloids pulled through a solution of A-DNA 
(used here as a monodisperse model polymer) with an optical tweezer. The experiments show a violation of the 
Stokes-Einstein relation based on the independently measured viscosity of the DNA solution: the drag force 
is larger than expected. We attribute this to the accumulation of DNA infront of the colloid and the reduced 
DNA density behind the colloid. This hypothesis is corroborated by a simple drift-diffusion model for the DNA 
molecules, which reproduces the experimental data surprisingly well, as well as by corresponding Brownian 
dynamics simulations. 



I. INTRODUCTION 

Complex fluids in general and colloid-polymer mixtures in 
particular are an ideal model system for studying of the struc- 
ture and phase behaviour of multicomponent systems. But 
they also play a large role in many technological processes 
such as oil recovery, food science 111, as well as in most bi- 
ological systems. For these systems, in addition to the equi- 
librium properties the dynamics is important. While until re- 
cently, research largely focused on bulk rheological proper- 
ties, it is now understood that interface effects play a signif- 
icant role. The structure of fluids in general changes in the 
vicinity of interfaces. In equilibrium situations this is known 
to lead to fluid structure mediated interactions, i.e., the so- 
called solvation force or depletion interactions. Theoretically, 
the concept of depletion interactions in colloid-polymer solu- 
tions has been extended to non-equilibrium situations |0, 
While in equilibrium, these interactions are short ranged (with 
an interaction range on the order of the particle diameter), they 
become long ranged in non-equilibrium situations. The rea- 
son for this are the long ranged structural changes in the fluid. 
These structural changes, i.e., an enhanced polymer density 
infront of the colloid, and a reduced polymer density in its 
back, also lead to an enhanced friction Et]. In addition, 
the structural changes in the vicinity of walls couple back to 
the rheologicalproperties there, e.g., leading to slip boundary 
conditions |0, 0]. This, on the other hand, leads to a drag 
reduction as compared to the expected Stokes drag for the 
bulk fluid. This work was inspired by dynamic light scatter- 
ing experiments on coUoids in solutions of macro molecules 

iilllii. 
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In order to elucidate the dynamics of colloids in polymer 
solutions in more detail, we perform highly controlled exper- 
iments on single isolated colloids in A-DNA solutions using 
optical tweezers. This technique allows us to move a col- 
loid through a highly mono-disperse polymer solution at a 
given velocity and to measure the drag force on the colloid 
with pN-resolution at the same time. For high DNA con- 
centrations we find a significantly higher drag force than pre- 
dicted by the Stokes-Einstein theory for the homogeneous so- 
lution, however, the force scales approximately linearly with 
the velocity. We compare the experimental results to a simpli- 
fied dynamic density functional theory (DDFT) for the non- 
interacting polymers in a flowing solvent and to Brownian dy- 
namics (BD) simulations. 



II. EXPERIMENT 

A. Materials and methods 

Fig.[T|illustrates the experimental set up, in which one col- 
loid is held by an optical trap surrounded by a polymer so- 
lution of A-DNA (obtained by New England BioLabs, Ger- 
many). An inverted microscope (Axiovert S 100 TV, Carl 
Zeiss, Jena, Germany) is used and the optical trap is real- 
ized with a diode pumped Nd:-YAG laser (1064 nm, 1 Watt, 
LCS-DTL 322; Laser 2000, Wessling, Germany). Its power 
is stabilized to achieve long-term stability. Additionally, the 
profile of the laser-beam was monitored. After passing an op- 
tical isolator, a quarter-wave plate is used in order to produce 
circularly polarized light to exclude effects due to reflection 
differences of the mirrors between the p- and s-polarisation 
of the laser light. The beam is expanded and coupled into 
the back aperture of the microscope objective (Plan-Neofluor 
100x1.30 Oil, Carl Zeiss, Jena, Germany). Video imaging 
and the optical position detection are accomplished by a dig- 
ital camera (KPF 120, Hitachi, Diisseldorf, Germany). The 
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FIG. 1: Illustration of the experimental setup. The colloid (grey) 
sourrounded by coils of DNA (red) is hold in an optical trap estab- 
lished by a photonic potential (magenta). 



optical stage is positioned in three dimensions with nanome- 
ter resolution using piezoactuators (P-5173CD, Physik Instru- 
mente, Karlsruhe, Germany). The sample cell consists of a 
closed chamber that can be flushed by a syringe pump with 
varying solutions. The viscosities are measured with an Ost- 
wald viscosimeter at the same temperature (25°C room tem- 
perature) and with the same A-DNA solutions as used in the 
experiment. 



B. Data analysis and calibration 

The position of the bead in the optical trap is determined by 
image analysis (Fig.|2^) llTll[l2ll . For that a sequence of im- 
ages is recorded (repetition rate 30 Hz) and analysed based on 
the Levenberg-Marquardt algorithm (Fig.|2t)). As fit function 

I = Io + A{1- ad) cxp{-d^) (1) 

is used, where {xq, yo) is the center position, r the optical ra- 
dius, A the amplitude of the profile relative to the background 
image intensity /o and a a constant to consider the dark dif- 
fuse ring around the colloid. A variation of r, A or a in an 
image sequence indicates the motion in z-direction. The cali- 
bration of the optical trap is based on Stokes' law for the pure 
solvent (here water) as described in detail elsewhere ||13(]. A 
typical force constant of the trap is 0.085 pN/nm correspond- 
ing to forces in the range between 0-50 pN, which can be 
determined with an accuracy of ±0.15 pN. 

In the experiments, no sign of irreversible adsorption of A- 
DNA molecules to the colloid was detected. While switching 
between different flow velocities, the same forces at the same 
speed were obtained in the range of our uncertainties and after 
flushing back to pure water the Stokes force of the colloid 
could be reproduced. 



FIG. 2: (a) microscope image of a single colloid (1.12/im radius) 
in the optical trap, (b) fit of the intensity distribution for the image 
shown in (a) using the Levenberg-Marquardt algorithm. 
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FIG. 3: Drag force F on colloids of radii Rc = 1.5 /im (boxes), 
1.12 /^m (circles), and 0.66 /im (diamonds) in pdna = 30 /xg/ml 
(full symbols) and 20 ^ig/ml (half filled symbols), as well as in pure 
water (open symbols) as a function of the pulling speed u. Also 
shown are linear fits to the data from which one can extract a density 
dependent viscosity via the Stokes-Equation However, as shown 
in Fig. [4] these viscosities are inconsistent as they depend on the 
colloid radius. 



C. Experimental results 

Fig. [3] shows the mesured force on colloids of sizes Rc = 
1.5 i-im, 1.12/im, and 0.66 /im as a function of the dragging 
velocity u in pure water as well as for DNA concentrations 
Pdna ~ 20 /ig/ml and pdna = 30 /ig/ml. Linear fits to the 
data analyzed in terms of the Stokes drag formula for a colloid 
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Pdna = 


O^g/ml 


20 ^ig/ml 


30 ^tg/ml 


Rc = 1.50 /im 


1.03 mPas 


1.36 mPas 


2.07 mPas 


Rc = 1.12 


1.00 mPas 


1.12 mPas 


1.30 mPas 


Rc ~ 0.66 /im 


1.01 mPas 


1.21 mPas 


1.37 mPas 


viscosimeter 


1.00 mPas 


1.09mPa.s 


1.14 mPas 



TABLE I: Viscosity as a function of DNA concentration extracted 
from linear fits to the data shown in Fig. [3] using the Stokes formula 
Eq. ([3) compared to the viscosities measured in a viscosimeter. From 
the viscosimeter data one obtains an intrinsic viscosity of the DNA 
of M = 1.55. See also Fig. H 

2.5 n < 1 < < n 

Tl = Tl„ (1 + [Tl] PoNA V/M) with [Tl] = 1.55 

* viscosimeter 
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FIG. 4: Viscosity as a function of DNA concentration extracted from 
linear fits to the data shown in Fig.|3]using the Stokes formula Eq. l|3) 
compared to the viscosities measured in a viscosimeter, c.f., Table|ll 
From the viscosimeter data one obtains an intrinsic viscosity of the 
DNA of [77] = 1.55. 

of diameter Rc in a DNA solution of viscosity t] 

^6nriRcu (3) 

yield the DNA concentration dependent viscosities in Table U 
shown in Fig. |4l The viscosities obtained in this manner 
are significantly larger than the viscosities measured in a vis- 
cosimeter for the same DNA solution and they increase with 
the colloid radius Rc- This is a strong indication that on top 
of the Stokes drag a second mechanism plays a role. 

Assuming the viscosities measured in the viscosimeter we 
can also extract an effective hydrodynamic radius i?eif from 
a fit of Eq. (O to the data in Fig. |3] The resulting i?eff as 
a function of pdna normalized to the nominal colloid radius 
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FIG. 5: Effective colloid radius i?eff extracted from the linear fits to 
the data in Fig.[3]relative to the nominal radius Rc as a function of 
the DNA density pdna- 



PDNA = 


A*g/ml 


20 pg/ml 


30 pg/ml 


Rc ~ 1.50 pm 


1.53 pm 


1.87 pm 


2.72 pm 


Rc = 1.12 Aim 


1.12 pm 


1.14pm 


1.27 pm 


Rc ~ 0.66 ^m 


0.66 pm 


0.73 pm 


0.79 pm 



TABLE II: Effective colloid radius R^ff extracted from the linear fits 
to the data in Fig. [3] 
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FIG. 6: Drag force F on colloids of radii Rc = 1.5 pm (boxes), 
1.12 pm (circles), and 0.66 pm (diamonds) in pdna = 30 pg/1 (full 
symbols) and 20 pg/1 (half filled symbols), as well as in pure water 
(open symbols) as a function of the pulling speed it. Also shown is 
the Stokes force Fs on the colloids as expected for the viscosities 
measured for the given DNA concentration in a viscosimeter, c.f., 
TableUandFig.E] 

Rc is shown in Table HIl and Fig.|5] Res/Rc increases with 
Pdna but it also varies with Rc'- the value for Rc = 1.5 pm 
is much larger than the values for the smaller colloids. If the 
increase of i?eff was due to adsorption of DNA molecules to 
the colloid we would expect the absolute increase of i?eff to a 
first approximation to be independent of Rc- 

The molecular weight of a A-DNA molecule is M ~ 
31.5 X lO^amu = 5.23 x 10"^^ pg and its contour length 
is about 16 pm, leading to a radius of gyration of roughly 
Rg = 0.5 pm. This corresponds roughly to a hydrodynamic 
radius of Rh ~ 0.662 i?c, = 0.33 pm. Based on the hydro- 
dynamic radius and the radius of gyration we get molecular 
volumes Vr = 1.5 x IQ'^'^ m\ and Vg = 5.2 x lO-^^mi^ 
respectively. For the highest DNA concentrations used in the 
experiment, i.e., pdna = 50 ptg/ml this corresponds to vol- 
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FIG. 7: The measured drag force J" on a colloid of radius 1.12 fim 
normalized to the velocity u as a function of the DNA concentration 
(open symbols). For the clarity of the presentation only a subset of 
the experimental data is shown. The data collapse demonstrates that 
F is proportional to it, however, the dependence on pdna is non- 
linear. The drag is also significantly larger than expected from the 
increased viscosity as measured in a viscosimeter (*). In pure water 
we obtain F/u — 0.021 pN/(/im/s) (dashed line). Also shown are 
simulation results for polymers with modified mobility (as explained 
in the text — full symbols). A fit for concentrations between and 
20/xg/mI (full line) highlights the nonlinearity in the density. 



ume packing fractions of 14% and 50%, respectively. With 
this, we can extract the intrinsic viscosity of the DNA solu- 
tion from the viscosimeter data and get — 1.55. Assuming 
the intrinsic viscosity for hard sphere suspensions [77] —5/2 
as predicted in [lAj, |15|] we get an effective hydrodynamic ra- 
dius of i?,, = 0.28 fim. 

In Fig.|6]we compare the measured drag forces on the col- 
loids with the predictions given by the Stokes equation 
assuming the viscosity values measured independently in the 
viscosimeter The data is normalized to the measurement in 
pure water. Apparently, in the DNA solution the drag force 
is much larger than expected for the increase in viscosity. In 
addition, the difference of the measured force to the predic- 
tion assuming Stokes' equation (|3]l with the viscosity values 
measured in the viscosimeter is larger for larger colloids. This 
is also apparent in Figs. |4] and |5] where the spurious viscosity 
and the effective hydrodynamic radius both increase not only 
with the DNA concentration but also with the colloid radius. 

Fig. I2] shows the drag force normalized to the pulling ve- 
locity w on a colloid of diameter Rc = 1.12 /_tm measured 
as a function of the DNA-density pdna- The data point at 
Pdna = is normalized to the Stokes drag force with the vis- 
cosity of water r^uj = 10^'^Ns/(m^). Clearly, the additional 
drag due to the presence of the DNA in the solution is not hn- 
ear in pdna and larger as the drag expected from the increased 
viscosity. The the errorbars for small velocities in Fig. Q are 
large due to the error in the velocity. 

In Fig.[8]the measured drag force on the colloid is compared 
to the Stokes drag force Fg in the DNA solution, see Eq. (O. 
The forces are to a good approximation linear in u. The differ- 
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FIG. 8: Drag force F on a colloid of radius 1.12 /im measured in 
pure water (Q) and in a DNA solution (□) as a function of the ve- 
locity It. The DNA concentration in (a) is 20 fig/ml and in (b) it is 
40 /ig/ml. The data is compared to the Stokes friction Fs in pure 
water (full line) and in the DNA solution Fg (dotted line) calculated 
from the measured viscosity, as well as to Fg plus the contribution 
Fp from the DNA jam in front of the colloid for Rp = 0.5 fim 
(dashed line) and for Rp —0.7 fim (dashed-doted line). Also shown 
is a fit to BD simulation results between O/xm/s and 50/^m/s (dash- 
dot-dotted line). 



ence between the drag forces for 20 /ig/ml DNA and for pure 
water is only significant for velocities larger than 400 /im/s. 
For 40 /ig/ml DNA we can observe a significant difference 
already at velocities larger than 100 /im/s. The experimen- 
tal resolution limits the measurment of the DNA induced drag 
force to relatively large velocities and to large DNA concen- 
trations. However, Fig.[8jb) clearly shows that the measured 
drag forces cannot be explained simply by an increased vis- 
cosity ?7dna of the solution. 



III. DRIFT-DIFFUSION MODEL 

That the drag force on the colloid cannot be explained sim- 
ply by the increased viscosity of the DNA solution and the 
dynamics of the DNA in the vicinity of the moving colloid 



5 



^ ^ ^ ^ ^ ^ 

^ ^ ^ ^ -5 

■U. ^ ^ l ^.i 

^ ^ ^ ^ / i 

^ ^ ^ i 

^ I I ^' i 

^ ^ ^ ^ ^ 

^ -a (I * 

Jj. ^ J> t J> 

J]- H. A 

^ J> ^ 4 ft 

-J ^1. 5. «?\iv & 

Ji. 4> 5. fi- V ' 

^ ^ ^ h \ h >. 

I ^ Si, h i> ^ 

^ ^ S^ h i, 

^ ji- 5. 5V 5V 



4 ft 




i> 
^ 

. ^ 

*S JJ. ^ ^ 

4). ^ 
^{ * * 4 



^ ^ ^ 
* ^ ^ 



^ -5. ^ 

■ ^ ^ ^ ^ 

4 {!. ^ ^ 

.4-'^ 4 4 4 S. ^ 

i 4 4 ^ 4 



FIG. 9: Flow field v(r) (arrows) around a moving colloid (full circle) 
with radius Rc in the frame of reference comoving with the colloid. 
Particles (e.g., DNA, open circle) with radius Rg can approach the 
colloid only up to a distance Rc + Rg (dashed circle). The flow field 
has a component normal to this circle which leads to an accumulation 
of particles infront and a depletion of particles in the back of the 
colloid. 



has to be taken into account, which is gouverned by the in- 
terplay of direct intermolecular interactions, hydrodynamics, 
and the internal degrees of freedom of a polymer chain. It 
has been shown, that an additional drag force due to the rear- 
rangement of solute particles in the vicinity of a dragged col- 
loid can be already obtained in a simple drift diffusion (DD) 
model lUt]. Here we employ the same model to calculate this 
additional drag force Fp and compare it to the experimental 
values. Within this model, the DNA molecules are idealized 
as mutually non-interacting particles with a finite hard core 
interaction radius with the colloid. The origin of this addi- 
tional force is illustrated in Fig. |9] Due to the repulsion of 
the DNA coils and the colloid, the center of mass of the DNA 
coils can approach the colloids surface only up to a distance of 
roughly Rg, i.e., the radius of gyration, which creates a forbid- 
den zone for the DNA. The solvent molecules, however, are 
much smaller and enter this zone, such that the solvent flow 
field has a component normal to the surface of the forbidden 
zone pointing inwards infront of the colloid and outwards be- 
hind the colloid. DNA coils advected with the solvent will 
therefore accumulate infront of the colloid and their density 
will be reduced in its back. This inhomogeneous DNA dis- 
tribution will lead to an inhomogeneous osmotic pressure and 
therefore to a force on the colloid. 

We will quantify the rearrangement of DNA molecules 
by calculating the average concentration C(r) (in units of 
molecules per volume) near the colloid, which, in the sim- 
ple drift-diffusion (DD) model described in |4i,|5il|, is given by 
the stationary solution of the Smoluchowski equation 



w-VC ^ DAC 



(4) 



in a frame of reference comoving with the colloid, with the 
solvent velocity field v and the DNA (zero concentration) dif- 
fusion constant D. In contrast to Ref. v(r) is not uniform 



but it is the solution of Stokes' equation for a sphere of ra- 
dius Rc translating with velocity u through a resting solvent 
(v — u for |r| — > oo) (see, e.g., Ref. lll6ll for details). 



v(r) + u 

3Rc 
4r 



3^2 



u) 



, (5) 



with the unit vector f = r/|r|. As we model the interac- 
tion potential between the colloid and the polymers as a hard- 
sphere potential, the centers of the polymers are excluded 
from a sphere of radius Rc + Rg around the center of the col- 
loid, see Fig.|9l Therefore, the boundary condition for Eq. (HI 
on this sphere is 



i^r ■ i)l^R^+R^ =0, 



(6) 



i.e., the DNA current j ^ vC — D VC normal to the sphere's 
surface has to vanish. Far from the colloid, the DNA density 
should be constant, i.e., C(r) — *■ Cq for |r| — > oo. As the mu- 
tual interaction between the DNA molecules is neglected, the 
local pressure on the colloid surface can be calculated from 
the ideal gas law p = ksTC. Integrating the local pressure 
over the surface yields the force Fp on the colloid. 

The solution of Eq. (HI depends only on the dimensionless 



velocity u* 



,Rc+Rg_ 
D 



(the Peclet number) of the colloid and 



the size ratio of the involved particles R* = Rc /{Rc + Rg)- 
For Rg/Rc 0, i.e., for R* 1, the surface of the for- 
bidden zone coincides with the colloid surface and the DNA 
behaves like a solvent molecule. In this limit the solution of 
Eq. dUi is C(r) = Co, the DNA molecules do not accumulate 
infront of the colloid and therefore the additional drag force 
Fp vanishes. 

With a hydrodynamic radius of Rh = 0.33 /im the Stokes- 
Einstein relation leads to a diffusion constant Z? w 6 x 
10^ ^'^ m^/s. Therefore, the smallest velocities in the experi- 
ments {u = 50 fim/s) correspond to Peclet numbers u* larger 
than 100 (depending on the colloid radius). While for small 
u* Eq. (HJi can be solved analytically in linear order in the 
Peclet number, we have to use numerical methods for the 
experimental velocities. To this end, we expand the density 
field C(r) in spherical harmonics up to order N and obtain a 
system of + 1 ordinary differential equations for the |r|- 
dependent expansion coefficients which we solved numeri- 
cally with AUTO 2000 iH (for details, see Ref. For 
large u* a fine numerical discretization (large N) is needed 
since the thickness of the region in front of the colloid in 
which the DNA density is enhanced decreases with u* while 
the density in this region increases with u*. We therefore use 
N = 100, which allows us to calculate reliable solutions up 
to u* w 100. For such high velocities the drag force Fp is 
well approximated by an affine function (see Fig.[TO]and also 

Ref. m 

Fp{R\u*)/[{Rc + Rgf ksTCo] = + u*. 

(7) 

which we use to extrapolate to velocities larger than u* = 
100, see Fig. [TOl The coefficients for fits to the numerical 
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FIG. 10: Numerically calculated drag force for the radii used in the 
experiment as a function of the Peclet number u* (symbols) and the 
affine fits to the data for 80 < u* < 100. 



Rc [fJ,m] 


p* Rc 

^ Rc+Rp 




mi 


0.66 


0.57 


3.31 


0.126 


1.12 


0.69 


2.43 


0.063 


1.50 


0.75 


1.76 


0.036 



TABLE III: Fitting coefficients for the extrapolation of the force 
Fp to large velocities (Peclet numbers) u* for the three colloid sizes 
used in the experiments. The radius of gyration is assumed to be 
Rg = 0.5 /xm. 

data in the range 80 < it* < 100 for the experimentally rele- 
vant values Rc = 1.50, 1.12, 0.66 /.jm and Rg = 0.5 /.(m are 
given in Table Hill In order to calculate Fp also for other val- 
ues of R* and to test the sensitivity of our result to variation 
of Rg we use the linear interpolation of the coefficients a{R*) 
and P{R*) 

a{x) = 8.146 - 8.427 a;, (8a) 
/3(x) = 0.411 - 0.502 a;. (8b) 

For the ideal gas law which we use to calculate the local 
osmotic pressure on the colloid surface we need the number 
density Co rather than the mass density pdna as used in the ex- 
periments. With the molecular weight M = 5.23 x 10^^^ /ig 
we get Co [m-3] = 1.91 x lO^^ pdna [/ig/ml]. 

Fig. |8] compares the forces on a colloid of size Rc = 
1.12 fj,m predicted by the DD model to the experimental data. 
Although the difference of the force at a concentration of 
Pdna = 20 /ig/ml in pure water (i.e., only the Stokes force) is 
hardly significant and within the experimental error explained 
by the enhanced viscosity of the DNA solution (Fig. |8ja)), 
theory and experiment are in better correspondence (in partic- 
ular for larger velocities) if one takes into account the force 
Fp for a DNA radius of Rg = 0.5 /im. For higher DNA con- 
centrations (e.g., pdna = 40/xg/ml in Fig. |8jb)), where the 
measured forces cannot be explained by the increased viscos- 
ity, Fp calculated from the DD model for Rg ~ 0.5 /im is 



FIG. 11: A cut through a part of the simulated system (28 /im x 10 
/im X 1.4 /tm. Polymers are shown in red and the colloid in blue. 
The arrow indicates the direction of motion of the colloid. Even 
though the density map and the density profiles shown in Fig. [12] 
clearly show an accumulation of polymers in front of the colloid and 
a depletion region behind it, this cannot be seen from a single snap- 
shot of the system. 



too small. A reasonable fit to the data can be obtained for 
Rg = 0.7 /im, however, for this large value the forces pre- 
dicted for the lower DNA concentrations (see, e.g., Fig.[8ta)) 
are too large. 

The experimentally measured drag force as a function of the 
DNA concentration as shown in Fig.|2]clearly shows a nonlin- 
earity indicating that interactions among the DNA molecules 
are relevant for concentrations larger than /jqna = 20 /ig/ml. 
The DD model presented in this section does not take these 
into account. Therefore, Brownian dynamics simulations are 
performed as presented in the following section. 

IV. BROWNIAN DYNAMICS SIMULATIONS 

We simulate the experiments with colloids of radius Rc ~ 
1.12 /im using a Brownian Dynamics (BD) method - also 
called Langevin Dynamics |17il . The polymers and the colloid 
are modeled as hard spheres with their respective radii. For the 
polymers we use the radius of gyration Rg = 0.5 /tm. As the 
polymers drag along most of the water contained in their vol- 
ume, the polymers are assigned the mass of water contained 
in a sphere with radius Rg. The number of polymers in the 
system is chosen such that the number density (rather than 
the mass density) of the experiments is reproduced. We use a 
rectangular simulation volume of 100 x 20 x 20 /im^ with 
periodic boundary conditions in all three directions. For a 
polymer concentration of 40 /ig/ml this corresponds to 30592 
polymers. 

The colloid is trapped in a moving parabolic potential 
V{r) = \ar^, mimicking the optical tweezer The poten- 
tial has a spring constant of a = 7.5 x lO^^pN/nm, which 
gives a better signal to noise ratio than the experimental value 
of 8.5 X lO^^pN/nm. Figure [TT] shows a snapshot of our sim- 
ulation setup. 

In BD, two most important aspects of hydrodynamics felt 
by the suspended particles are taken into account, namely the 
Stokes friction and the Brownian motion. Correspondingly, 
this is done by adding to a molecular dynamics simulation 
two additional forces. The Langevin equation describes the 
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motion a Brownian particle with radius R at position r(t) as 



ir{t) = 6 7r?7i?r(i) +Frand(t) +Fext(r,t), 



(9) 



where the first term models the Stokes friction in a solvent 
of viscosity ry, Fext(r, t) is the sum of all external forces like 
gravity, forces excerted by other suspended particles, and, for 
the colloid, the optical trap. Frand(i) describes the thermal 
noise which gives rise to the Brownian motion. The random 
force on different particles is assumed to be uncorrected, as 
well as the force on the same particle at different times. It is 
further assumed to be Gaussian with zero mean. The mean 
square deviation of the Gaussian (i.e., the amplitude of the 
correlator) is given by the fluctuation-dissipation theorem as 



(iFrandP) = nnrjRkBT. 



(10) 



In order to reduce computation time, physical quantities are 
rescaled: the simulation is carried out at a lower temperature. 
To compensate for this, the viscosity of the fluid as well as all 
energies are scaled by the same factor In the present simula- 
tions, temperature is scaled down by a factor of 37500. This 
scaling leaves the diffusion constant as well as the relative im- 
portance of diffusion and motion caused by external forces 
unchanged. However, it allows for a much larger time step (in 
this case 60/is) [Ts^. 

From the simulation data, it is possible to measure the ef- 
fective polymer concentration around the dragged colloid. To 
accomplish this, about 600 snapshots of the simulation are 
taken: Each snapshot is moved, so that the position of the 
colloid coincides, and the probability for a certain space to 
be occupied by a polymer is calculated by averaging over all 
snapshots (200 x 100 bins are used). The picture in Fig. [T2h 
depicts a typical example of a density map for a concentra- 
tion of 20 ng/ml and u ~ 50 fim/s. It can be observed that 
there is a region of high polymer concentration infront and 
at the sides of the colloid whereas the region behind the col- 
loid is not yet completely refilled by the polymers. Figure[T2b 
shows the normalized average density of polymers in the di- 
rection of motion for concentrations of 5/^tg/ml and 20 /ig/ml 
and v=50 fim/s. Infront of the colloid, a sharp peak can be ob- 
served. In addition, on observes density oscillations which are 
typical for hard sphere suspensions. For high polymer concen- 
trations, the probability is close to 1 to find a polymer infront 
of the colloid. The region behind the colloid is almost clear of 
polymers, because polymers did not yet have enough time to 
relax into this region again. 

BD is widely used to simulate suspensions (e.g. |[l9i l20l 
I2TI1 ) because it is well understood, not difficult to implement, 
and needs much less computational resources than a full sim- 
ulation of the fluid. However, this simulation method does 
not resolve more complicated hydrodynamic phenomena, in 
particular long-ranged hydrodynamic interactions among the 
polymer particles and in particular between the polymers and 
the colloid. As a result both, the friction and the thermal force 
excerted on a particle are independent of the position and ve- 
locity of the other particles. In the case of a polymer trapped 
infront of the moving colloid (i.e., if it is dragged along by the 
colloid), this leads to a considerable overestimation of the fric- 
tion felt by the colloid: in a BD simulation it feels the friction 




FIG. 12: (a) Polymer density around the colloid averaged over 
600 snapshots of the system at a concentration of 20 ^g/ml and 
u — 50 /im/s. Lighter colors denote higher polymer densities, (b) 
Normalized average polymer density on a line through the colloid 
center in the direction of motion for concentrations of 5 /ig/ml (cor- 
cles) and 20 /ig/ml (boxes) and velocity u = 50 /im/s. Polymers 
accumulate infront of the colloid and the concentration in the back 
is reduced due to the finite Peclet number of the polymers. Lines 
connecting datapoints are guides to the eye. Also visible are density 
oscillations infront of the colloid, which are characteristic for hard 
sphere systems. The vertical lines indicate the size of the colloid 
(±1.12 /im). 



on itself and the dragged polymer, whereas in reality the fric- 
tion on the colloid would be significantly reduced because it is 
in the polymer's slipstream. In addition, as argued in Ref. iH, 
the hydrodynamic interaction of the colloid with the polymers 
is of importance. The solvent molecules flow around the col- 
loid and the resulting flow field advects the polymers around 
the colloid (see Eq. (|5]l), therefore significantly reducing the 
polymer jam infront of the colloid. To compensate for both of 
these effects we reduced the friction felt by the polymers, i.e., 
increase their diffusion constant. This does not significantly 
influence the density distribution outside the interaction re- 
gion where the polymers are homogeneously distributed. 

Our model reproduces very well the linear relation between 
drag force and the drag velocity for different polymer con- 
centrations. As higher drag velocities require a larger system 
(even with periodic boundary conditions) and short numerical 
timesteps, we are limited to about fifty micrometers per sec- 
ond by the available computational resources and time. How- 
ever, this is not a problem because of the linearity of the drag 
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FIG. 13: Effective viscosity (calculated with Eq. ^ from drag force 
measured in the simulations) versus drag velocity for different poly- 
mer concentrations. For velocities larger than about 10 fim/s, the 
effective viscosity is independent of the velocity as expected. For 
low velocities, the drag force increases. The horizontal line indicates 
the viscosity of the solvent water. 



force with respect to velocity for higher velocities as observed 
in the experiments. Fig.[8]includes a linear fit to the simulation 
data for velocities between 0/im/s and 50/^im/s, which agrees 
well with the experimental data. 

As in the experiments, with Eq. (O the measured drag 
force can be interpreted in terms of an effective viscosity. 
Fig. [13] shows the dependency of this viscosity on the veloc- 
ity for pure water and for polymer concentrations of 20^g/ml 
and 40/ig/ml. Even for a polymer concentration of 40/^(g/ml 
the viscosity stays constant for drag velocities greater than 
30/im/s, but for very low velocities, in all three cases we 
find that the drag force decreases less then linearly. For 
such low drag velocities the colloid is very near the trap cen- 
ter and thermal fluctuations about the average position are 
large. An energy of /cbT corresponds to a displacement of 
0.3 ^m, whereas the displacement calculated from the drag 
force in water at u = 1 /^jm/s corresponds to a displacement of 
X = 0.28^m. However, these fluctuations should be roughly 
symmetric with respect to the average colloid position within 
the trap. However, the force calculated numerically within 
the framework of the DD model shown in Fig. [TO] analyzed in 
terms of an effective viscosity would also give an increased 
viscosity for smaller velocities: the slope is larger for smaller 
values of u*. Therefore we conclude, that the apparent in- 
crease in viscosity as shown in Fig. [13] is at least partially 
related to the formation of the polymer jam infront of the 
colloid. However, the dataset for pure water (i.e., without 
polymers) indicates, that at very low velocities, fluctuations 
do play a role. 

For a fixed drag velocity and varying polymer concentra- 
tions, the linear relation between force and concentrations 
for low polymer concentrations and the stronger nonlinear 
increase for higher concentrations can also be observed as 
shown in Fig. [7] The increase in drag force is overestimated 
as expected if the mobility of the polymers is not increased 
(as described above). When the mobility of the polymers is 
increased by a factor of 14.3, however, quantitative agreement 
is achieved for concentrations up to 40 /ig/ml. For even higher 
concentrations, the drag force increases stronger than in the 



experiments. This is because at high concentrations, a hard 
sphere suspension has a stronger tendency to jam than a real 
polymer suspension: polymers are deformable such that, in 
contrast to hard spheres, they can squeeze past each other The 
simulation results for increased polymer mobility indicate that 
the friction felt by a pair of a colloid and a dragged polymer 
in the experiment is far less than the sum of the friction forces 
on a single polymer and colloid, respectively. 



V. DISCUSSION 

The experiments clearly show that the drag force on col- 
loids pulled through a solution of A-DNA with an optical 
tweezer cannot be explained by the Stokes force for the vis- 
cosity of the solution. It is much higher and increases non- 
linearly with the DNA-concentration but approximately lin- 
early with the drag velocity. The resolution of the force mea- 
surment is ±0.15 pN and therefore the difference between the 
forces on the colloid in the solution and in pure water can 
be measured only for large concentrations or, in the case of 
smaller concentrations only for large velocities. This limits 
the overlap between the DD theory and the experiments to the 
lowest concentration of pdna = 20 /ig/ml used in the experi- 
ment. 

At this concentration, the volume fraction of the solution 
taken by the DNA coils (based on Rg = 0.5 /im) is 0.1 and 
DNA-DNA interactions can be neglected. The agreement be- 
tween theory and experiment is very good. For the highest 
concentration used in the experiments pdna = 50 /ig/ml the 
volume fraction is 0.5 and DNA-DNA interactions have to be 
taken into account. 

In the BD simulations DNA-DNA interactions have been 
taken into account. As in the experiments and in the DD the- 
ory, the drag force increases linearly with the velocity, at least 
for large velocities. The nonlinear increase of the drag force 
with the DNA concentration is also observed in the BD sim- 
ulations. However, since the BD simulations could not take 
into account hydrodynamic interactions between the colloid 
and the DNA molecules, the drag force is significantly over- 
estimated. The hydrodynamic flow field around the colloid as 
shown in Fig. [9] has a component normal to the direction of 
motion of the colloid which efficiently reduces the number of 
DNA molecules accumulated infront of the colloid, while it 
fills the depletion zone at its back. As demonstrated by the 
excellent agreement with the experimental data, this hydrody- 
namic interaction between the colloid and the DNA has the 
same effect as increasing the DNA mobility, here by a factor 
of 14.3. This is consistent with results of BD simulations of 
soft particles in a background flow, where the drag force was 
reduced by a factor of 10 by this hydrodynamic effect js]]. 

The agreement between theory and experiment for low con- 
centrations but high velocities is particularly remarkable con- 
sidering the simplicity of the models which neglect many as- 
pects of hydrodynamic interactions or include them in a rather 
simplistic way (e.g., in terms of an increased viscosity in the 
DD model). The internal structure of the DNA molecules is 
neglected in both, the DD model as well as the BD Simula- 
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tions. Apparently, the deformation of the DNA coils in the 
accumulation region infront of the colloid is not significant. 
This can be only explained by the proximity of the stagnation 
point for the solvent flow which limits the normal component 
of the the solvent velocity field. At this point the shear rate 
is also smaller than at the colloid's side, where it goes up to 
over 700 1/s for u = 500 /^m/s. With a relaxation time on the 
order of 0. 1 s ll22ll this amount to a Weissenberg number Wi 
on the order of 70, i.e., much larger than one. Therefore the 
polymercoils should be significantly distorted and viscoelastic 
effects are to be expected. 

In summary, we present a first direct experimental obser- 
vation of jamming-induced drag-enhancement on colloids in 
polymer solutions. First theoretical approaches to this prob- 



lem as well as the Brownian dynamics simulations presented 
in this paper contain many simplifications which need to be 
addressed in the future. 
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